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LONG-TERM  GOAL 

Our  long-term  goal  is  to  contribute  to  the  development  and  evaluation  of  an  operational  evaporation 
duct  model  for  use  in  microwave  propagation  prediction. 

OBJECTIVES 

Our  first  objective  is  to  develop  a  metric,  involving  one  or  more  parameters,  to  replace  evaporation 
duct  height  as  a  fundamental  basis  for  comparison  and  evaluation  of  evaporation  duct  models.  Because 
electromagnetic  propagation  is  determined  by  vertical  gradients  of  modified  refractivity  (M)  and 
because  these  gradients  are  more  uniquely  determinant  of  propagation  conditions  that  duct  height  alone 
(Dockery  1988,  Reilly  and  Dockery  1990,  Gehman  2000),  it  is  felt  that  evaporation  duct  height  is  often 
inadequate  for  accurately  describing  propagation  conditions.  Our  second  objective  is  to  use 
experimental  data  to  evaluate  various  LKB-based  (where  LKB  refers  to  Liu  et  al  1979)  evaporation 
duct  models  (e.g.,  the  BYC  model  described  in  Babin  et  al  1997,  the  NRL-Monterey  Model, .etc.)  to 
determine  the  best  model(s)  using  the  objective  criterion  from  the  first  objective.  Our  third  objective  is 
to  evaluate  the  possible  effects  of  issues  raised  during  the  Naval  Postgraduate  School  Evaporation  Duct 
Workshops  on  the  LKB  approach  to  evaporation  duct  models. 

APPROACH 

Field  experiments  were  conducted  in  which  the  R/V  CHESSIE  was  deployed  off  the  coast  of  Wallops 
Island,  VA.  A  surface  buoy  was  developed  by  JHU/APL  (John  Rowland  and  Charles  Etheridge)  for 
making  fine-scale  measurements  of  temperature  and  humidity  profiles  between  about  0.1  and  0.8  m 
above  the  sea  surface  using  a  sensor  package  that  travels  up  and  down  a  mast  with  a  one-way  time  of 
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seven  minutes.  This  mast  is  mounted  on  a  catamaran  to  reduce  wave  effects.  This  catamaran  is  at  the 
end  of  a  46  m  tether  attached  to  the  stem  and  upwind  side  of  the  R/V  CHESSEE.  The  measurements 
made  by  this  profiling  buoy  were  compared  with  profiles  generated  by  four  LKB-based  evaporation 
duct  models.  These  models  use  air  temperature,  humidity,  wind  speed,  pressure  and  sea  surface 
temperature  (SST)  as  input.  These  inputs  were  measured  by  sensor  packages  mounted  on  masts  on  the 
R/V  CHESSIE  (data  collected  on  the  upwind  mast  to  reduce  ship  contamination)  and  by  an  SST  buoy 
at  the  end  of  a  46  m  tether  attached  to  the  bow  and  upwind  side  of  the  R/V  CHESSIE.  This  SST  buoy 
was  about  15  m  from  the  profiling  buoy  when  deployed.  These  measurements  were  saved  as  five- 
minute  averages.  Because  the  profiling  buoy  individual  profiles  last  seven  minutes  and  the  model 
inputs  were  five-minute  averages,  the  minimum  time  for  contemporaneous  measurements  is  35- 
minutes.  In  order  to  reduce  the  effects  of  turbulence  in  the  data  and  because  the  LKB  models  were 
intended  to  produce  ensemble  average  profiles,  70-minute  averages  of  all  data  were  used. 

WORK  COMPLETED 

Software  for  reading,  analyzing  and  plotting  the  data  was  written.  Several  alternatives  to  evaporation 
duct  height  as  a  criterion  for  model  comparison  were  considered.  Three  of  these  criteria  were 
examined  in  detail  and  were  used  to  compare  model  generated  profiles  with  data  from  the  JHU/APL 
profiling  buoy:  1)  The  mean  M  slope  root-mean-square  (RMS)  difference  criterion  is  determined  by 
taking  the  mean  of  all  the  M  slopes  between  each  contiguous  pair  of  data  points  and  each  pair  of 
model-derived  points  over  the  height  range  of  profiling  buoy  data.  For  each  set  of  data,  there  is  then  a 
mean  data  M  slope  and,  for  each  model,  a  mean  model  M  slope.  The  difference  between  these  means 
for  each  set  of  profiles  is  then  calculated,  this  difference  is  squared,  and  the  square  root  of  the  mean  of 
all  these  squared  differences  for  all  data  sets  is  the  final  criterion.  2)  The  RMS  M  value  criterion  is 
determined  by  taking  the  M  at  each  height  in  a  profile,  subtracting  the  mean  M  of  that  profile  and 
taking  the  RMS  difference  between  these  M  values  for  a  model  and  the  data.  3)  The  RMS  M  slope 
difference  criterion  is  determined  by  making  a  log-linear  curve  fit  to  the  buoy  profile,  finding  the  M 
slopes  for  each  0.02  m  segment  of  this  curve-fit  profile  and  the  corresponding  altitude  segment  in  a 
model  profile,  and  finding  the  mean  RMS  difference  for  all  the  profiles. 

Four  LKB-based  evaporation  duct  models  were  compared.  They  include  the  NPS  model  (Frederickson 
et  al  2000),  the  NWA  model  (written  by  Kurt  Krai  of  the  Naval  Warfare  Assessment  Station  and  based 
on  Liu  and  Blanc,  1984),  the  NRL  model  (written  by  John  Cook  and  based  on  Liu  and  Blanc,  1984), 
and  the  BYC  model  (Babin  et  al  1996).  The  primary  difference  among  these  models  is  that  they  use 
different  'F  functions.  These  functions  are  used  in  Monin-Obukhov  similarity  theory  (Monin  and 
Obukhov  1954)  to  account  for  differences  in  temperature,  humidity  and  wind  profiles  from  those  in 
neutral  stability.  Recall  that  Monin-Obukhov  similarity  theory  is  considered  valid  within  the  surface 
layer  for  both  stable  and  well-mixed  conditions  (Stull  1994),  although  the  equations  were  derived 
primarily  from  experiments  over  land.  Fairall  et  al  (1996)  used  LKB  surface  layer  theory  and  other 
results  to  extend  Monin-Obukhov  similarity  to  conditions  over  water,  although  they  encountered 
mostly  unstable  conditions.  Therefore,  the  'F  functions  used  by  all  the  models  for  stable  conditions 
were  derived  from  experiments  over  land.  Furthermore,  wave  effects  are  not  explicitly  parameterized 
by  any  of  these  models. 

A  total  of  421  M  profiles  were  collected  during  field  tests  in  November  1998,  the  Y2K  Wallops 
Microwave  Propagation  Measurement  System  (MPMS)  Experiment  in  April-May  2000  and  other  field 


experiments  in  August  and  September  2000.  Examples  of  single  profiles  are  shown  in  Figs.  1  and  2, 
where  the  air-sea  virtual  temperature  differences  were  -4.7C  and  3.2C,  respectively.  Recall  that  it 
takes  7  min  for  the  sensor  to  transit  the  mast.  Note  how  there  is  more  scatter  for  unstable  than  for 
stable  conditions.  Also,  note  the  waviness  in  the  profiles.  While  we  did  not  have  wave  measurements 
for  corroboration,  it  is  possible  that  this  waviness  in  the  profiles  is  induced  by  waves  in  the  ocean.  Fig. 
3  shows  a  comparison  between  the  model  profiles  and  the  April  30  averaged  profiling  buoy  data.  Fig.  4 
is  similar  to  Fig.  3  but  for  the  September  22  unstable  case.  Fig.  5  shows  just  the  model  profiles  for  this 
case  to  illustrate  how  they  diverge  at  altitudes  above  1  m,  and  also  to  illustrates  why  a  model  duct 
height  is  better  determined  using  dM/dz=0  instead  of  the  minimum  M  value. 
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Figure  1  (Unstable)  Figure  2  (Stable) 


Figure  3  Figure  4 
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Figure  5 


RESULTS 


A  variety  of  atmospheric  conditions  were  observed,  with  air-sea  virtual  temperature  differences  ranging 
from  -4.7C  to  3.2C  and  wind  speeds  ranging  from  2.8  to  8.6  m/s.  Using  profiles  adjacent  in  time,  a 
total  of  thirty-four  70-minute  averages  were  obtained  that  had  contemporaneous  model  input  data  from 
R/V  CHESSIE  and  the  SST  buoy.  These  34  data  sets  were  examined  using  the  comparison  criteria 
above.  While  70-minute  averages  were  used  for  the  comparison,  the  overall  individual  profile  shapes 
were  fairly  repeatable. 


When  the  RMS  M  slope  difference  criterion  (number  3  above)  is  used,  the  NPS  model  performs  the 
best  for  stable  and  near-neutral  conditions,  while  the  BYC  model  performs  the  best  for  unstable 
conditions.  Overall,  the  RMS  M  slope  differences  were  on  the  order  of  2  M/m  for  stable  conditions, 
less  for  near-neutral  conditions  and  more  for  unstable  conditions.  A  plot  of  the  RMS  M  slope 
differences  for  the  NPS  model  versus  air-sea  virtual  temperature  difference  is  shown  in  Fig.  6.  Plots 
for  the  other  models  show  a  similar  variation  with  stability.  When  the  RMS  M  difference  after  mean 
removal  criterion  is  used  (number  2  above),  the  BYC  model  performs  the  best  for  all  conditions. 
Overall,  the  RMS  M  value  differences  were  on  the  order  of  2  M  units  for  unstable  conditions  and  less 
than  1  M  unit  for  stable  and  near-neutral  conditions.  A  plot  of  the  RMS  M  difference  for  the  NPS 
model  versus  air-sea  temperature  difference  is  shown  in  Fig.  7.  When  the  criterion  comparing  the 
mean  slope  of  the  entire  profile  between  the  model  and  the  data  (number  1  above)  was  used,  the  results 
were  the  same  as  those  for  the  RMS  M  slope  difference  (number  3  above).  The  mean  M  slope  RMS 
differences  were  on  the  order  of  3  M/m  for  unstable  conditions  and  1  M/m  for  near-neutral  and  stable 
conditions.  What  these  errors  mean  in  terms  of  propagation  is  a  question  yet  to  be  resolved,  primarily 
because  these  errors  are  in  the  region  below  1  m. 
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Figure  6 


Figure  7 


Because  there  are  cases  in  which  the  models  do  a  relatively  poor  job  of  matching  profile  shape,  it 
seems  likely  that  there  are  some  phenomena  occurring  in  the  region  below  1  m  that  is  missing  from  the 
models.  Observations  of  the  waviness  about  a  mean  profile  shape  appearing  in  both  individual  profiles 
and  in  the  70-minute  averages  lead  us  to  suspect  that  this  phenomena  is  primarily  the  effects  of  wave 
motion.  Unfortunately,  we  do  not  have  coincident  wave  measurements  to  corroborate  that  this  is  from 
wave  motion.  It  is  also  interesting  to  note  that,  while  the  model  profiles  often  look  similar  to  one 
another  below  1  m,  there  may  be  quite  a  divergence  at  higher  altitudes  (Fig.  5).  It  appears  that  the 
different  T  functions  used  by  the  models  have  an  important  effect  in  this  region.  It  would  therefore  be 


necessary  to  resolve  the  issue  of  which  ¥  function  is  best  for  stable  and  unstable  conditions  over  the 
ocean.  One  can  conclude  that  a  profiling  buoy  capable  of  measuring  higher  altitudes  and  measuring 
wave  effects  would  be  valuable  in  resolving  these  questions. 

It  has  been  observed  that  when  evaporation  duct  height  is  obtained  from  LKB-based  models,  the  output 
is  more  sensitive  to  input  under  stable  conditions  than  under  unstable  conditions.  Furthermore,  while 
the  original  functions  for  stable  and  unstable  conditions  were  developed  using  the  results  of 
experiments  over  land,  the  TOGA-COARE  experiments  resulted  in  modifications  to  the  'F  functions 
for  unstable  conditions.  Therefore,  one  might  expect  that  the  ¥  functions  would  perform  better  under 
unstable  conditions.  However,  we  have  the  interesting  result  that  any  one  of  the  four  LKB-based 
models  provides  better  results  for  stable  than  unstable  conditions.  The  difference  is  that  we  are 
comparing  profile  curvature  instead  of  duct  height.  Because  electromagnetic  propagation  models 
determine  propagation  using  M  slope,  this  should  be  a  more  reliable  comparison  than  duct  height. 

These  results  should  be  interpreted  with  caution  because  of  the  small  size  of  the  data  set.  While  the 
small  number  of  cases  makes  it  difficult  to  assign  a  clear  winner  model,  the  analytical  concept  used 
here  has  been  shown  to  be  a  valid  way  of  comparing  model  profile  curvatures.  The  small  altitude  range 
of  the  JHU/APL  profiling  buoy  also  limited  us.  This  altitude  range  of  about  0. 1-0.8  m  is  a  very  small 
fraction  of  the  0-40  m  range  produced  by  the  models.  A  new  buoy  is  being  developed  that  will  cover 
altitudes  up  to  about  8  m.  Applying  the  analytical  techniques  developed  in  this  study  to  an  expanded 
data  set  using  the  new  buoy  should  provide  a  more  robust  answer  to  the  question  of  evaporation  duct 
model  performance.  Obtaining  this  answer  becomes  even  more  important  as  radar  frequencies  go  from 
S  band  to  X  band. 

IMPACT/APPLICATION 

Operations  of  radar  and  communications  may  be  severely  affected  by  refractive  index  gradients  in  the 
lower  troposphere,  especially  over  the  ocean.  For  accurate  microwave  propagation  prediction,  there  is  a 
critical  need  for  an  evaporation  duct  model  that  provides  realistic  fine-scale  vertical  refractivity 
gradients  using  a  minimum  of  surface  layer  measurements.  Once  such  application  would  be  using  data 
from  the  MORIAH/SMOOS(R)  system  for  propagation  prediction.  Once  reliable  meteorological 
modeling  has  been  achieved,  simplifications  can  be  made  for  operational  purposes.  The  origin  and 
magnitude  of  any  errors  resulting  from  these  operational  simplifications  can  then  be  determined.  In 
addition,  the  resulting  evaporation  duct  model  may  eventually  be  used  with  the  near-surface  grid  point 
outputs  of  a  mesoscale  model  such  as  the  US  Navy’s  Coupled  Ocean/Atmospheric  Mesoscale 
Prediction  System  (CO AMPS)  to  provide  mesoscale  refractivity  forecasts. 

TRANSITIONS 

As  mentioned  under  Impacts  above,  the  results  of  this  effort  should  lead  to  useful  operational  products 
for  future  use  by  the  US  Navy.  For  example,  these  products  would  provide  a  way  of  using  data  from 
the  SMOOS (R)/MORIAH  system  for  propagation  prediction  of  evaporation  duct  effects. 

RELATED  PROJECTS 


None  at  this  time. 


REFERENCES 


Babin,  S.  M.,  G.  S.  Young,  and  J.  A.  Carton,  1997:  A  new  model  of  the  oceanic  evaporation  duct.  J. 
Appl.  Meteor.,  36, 193-204. 

Dockery,  G.  D.,  1988:  Modeling  electromagnetic  wave  propagation  in  the  troposphere  using  the 
parabolic  equation.  IEEE  Trans.  Antennas  Propag.,  36, 1464-1470. 

Fairall,  C.  W.,  E.  F.  Bradley,  D.  P.  Rogers,  J.  B.  Edson,  and  G.  S.  Young,  1996:  Bulk  parameterization 
of  air-sea  fluxes  for  Tropical  Ocean-Global  Atmosphere  Coupled-Ocean  Atmosphere  Response 
Experiment.  J.  Geophys.  Res.,  101, 3747-3764. 

Frederickson,  P.  A.,  K.  L.  Davidson  and  A.  K.  Goroch,  2000:  Operational  Bulk  Evaporation  Duct 
Model  for  MORIAH,  Naval  Postgraduate  School  Report,  Draft  Version  1.2,  May,  70  pp. 

Gehman,  J.  Z.,  2000:  Importance  of  Evaporation  Duct  Stability  in  Propagation-sensitive  Studies, 
JHU/APL  Technical  Report  A2A-00-U-3-008,  June,  8  pp. 

Liu,  W.  T.,  and  T.  V.  Blanc,  1984:  The  Liu,  Katsaros  and  Businger  (1979)  Bulk  Atmospheric  Flux 
Computational  Iteration  Program  in  FORTRAN  and  BASIC,  NRL  Memorandum  Report  5291,  8  May. 

Liu,  W.  T.,  K.  B.  Katsaros,  and  J.  A.  Businger,  1979:  Bulk  parameterization  of  air-sea  exchanges  of 
heat  and  water  vapor  including  the  molecular  constraints  at  the  interface.  J.  Atmos.  Sci.,  36, 1722- 
1735. 

Monin,  A.  S.,  and  A.  M.  Obukhov,  1954:  Basic  laws  of  turbulent  mixing  in  the  atmosphere  near  the 
ground.  Tr.  Akad.  Nauk.  SSSR  Geofiz.  Inst.,  24(151),  163-187. 

Reilly,  J.  P.,  and  G.  D.  Dockery,  1990:  Influence  of  evaporation  ducts  on  radar  sea  return.  IEE  Proc., 
PartF,  137(2),  80-88. 

Stull,  R.  B.,  1994:  An  Introduction  to  Boundary  Layer  Meteorology.  Kluwer  Academic  Pub.,  Norwell, 
MA,  666  pp. 

PUBLICATIONS 

Babin,  S.  M.  and  Dockery,  G.  D.,  2001:  Comparison  of  LKB-based  evaporation  duct  models  with  in 
situ  atmospheric  data.  Presented  at  the  URSI  National  Radio  Science  Meeting,  Boulder,  CO,  8 
January.  Preprints,  2001  USNC/URSI  National  Radio  Science  Meeting,  90. 

Babin,  S.  M.  and  Dockery,  G.  D.,  2001:  LKB-based  evaporation  duct  model  comparison  to  buoy  data. 
Submitted  to  Journal  of  Applied  Meteorology,  March  2001. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Pu^  reporting  ,0, »»  c^oono,  ^  ^ 

Papenworit  Radueaon  Preset (070W188)  Waahmgton.  OC  20501  .nnoPRC 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. - - - 1 ,  n ATFfi  COVERED  (From  -  To) 


1.  REPORT  DATE  (DD-MM-YYYY) 

30-03-2001 _ _ 

4.  TITLE  AND  SUBTITLE 


1 2.  REPORT  DATE 


Development  and  Experimental  Evaluation  of 
Oceanic  Evaporation  Duct  Models  Based  on 
the  LKB  Approach 


1 3.  DATES  COVERED  (From  -  To) 

Mav  1999  -  March  2001 

Sa.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

N000149910796 

5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Babin,  Steven  M. 
Dockery,  G.  Daniel 


|5d.  PROJECT  NUMBER 


|5e.  TASK  NUMBER 


|5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

The  Johns  Hopkins  University 
Applied  Physics  Laboratory 
11100  Johns  Hopkins  Road 

Laure 1  m  2b723-6099 _ 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 
Dr.  Scott  Sandgathe 
Office  of  Naval  Research 
Room  428,  800  N.  Quincy  Street 
Arlington,  VA  22217 

12.  DISTRIBUTION  AVAILABILITY  STATEMENT 

Unlimited 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

Final 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 


(See  attached) 


"  'Evaporation  duct,  electromagnetic  propagation,  marine  boundary  layer 


16.  SECURITY  CLASSIFICATION  OF: 
a  REPORT  lb.  ABSTRACT  |c.  THIS  PAGE 

u  U  u 


117  LIMITATION  OF  1 18.  NUMBER  1 19a.  NAME  OF  RESPONSIBLE  PERSON 

■  abstract  of  pages  Steven  M.  Babin _ 

UU  6  79tT"TELEP0NE7IuMBER"(//ic/lKJa"ariaeo5JJ"" 

240-228-6704 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI-Std  Z39-18 


Development  and  Experimental  Evaluation  of  Oceanic  Evaporation  Duct  Models  Based  on  the 
LKB  Approach 

Comparison  is  made  between  four  evaporation  duct  models  based  on  the  surface  layer  theory  of 
Liu,  Katsaros  and  Businger.  These  models  are  used  to  generate  modified  refractivity  (M)  profiles 
that  are  then  compared  with  those  determined  from  time-averaged  atmospheric  data  measured  by 
the  JHU/APL  profiling  buoy.  Model  inputs  are  derived  from  measurements  from  masts  on  the  R/V 
CHESSIE  and  a  tethered  sea-surface  temperature  (SST)  buoy.  Because  electromagnetic 
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employed  to  compare  the  curvature  of  the  model  profiles  with  that  of  the  profiles  measured  by  the 
profiling  buoy.  When  root-mean-square  M  slope  difference  between  model  and  a  curve  fit  to  the 
data  is  used  as  the  comparison  criterion,  the  NPS  model  performs  the  best  for  stable  and  near¬ 
neutral  conditions,  while  the  BYC  model  performs  the  best  for  unstable  conditions.  When  root- 
mean-square  M  difference  (after  mean  removal)  is  used  as  the  criterion,  then  the  BYC  model 
performs  the  best  for  all  conditions.  It  is  also  apparent  that  the  model-derived  profiles  may  be 
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